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Hunting for dark halo substructure using 
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ABSTRACT 

Dark halo substructure may reveal itself through secondary, small-scale gravitational 
lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we 
explore the prospects of using Very Long Baseline Interferometry (VLBI) observations 
of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions 
produced by various forms of dark substructures in the 10 3 -10 8 Mq mass range. We 
present lensing simulations relevant for the angular resolutions attainable with the 
existing European VLBI Network (EVN), the global VLBI array, and an upcoming 
observing mode in which the Atacama Large Millimeter Array (ALMA) is connected 
to the global VLBI array. While observations of this type would not be sensitive to 
standard cold dark matter subhalos, they can be used to detect more compact forms of 
halo substructure predicted in alternative structure formation scenarios. By mapping 
w 5 strongly lensed systems, it should be possible to detect or robustly rule out 
primordial black holes in the 10 3 -10 6 Mq mass range if they constitute > 1% of the 
dark matter in these lenses. Ultracompact minihalos are harder to detect using this 
technique, but 10 6 -10 8 Mq ultracompact minihalos could in principle be detected if 
they constitute > 10% of the dark matter. 
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1 INTRODUCTION 

A generic prediction of the standard cold dark matter 
(CDM) scenario is that a substantial fraction o f the total 
mass of galaxy-sized dark matt er halos (~ 10%; iGao et ail 
l201ll : iMacieiewski et all |201ll ) should be in the form of 
bound substructures (a.k.a. subhalos or subclumps) left over 
from the process of hierarchical assembly. The fact that the 
number of substructures seen in CDM simulations greatly 
outnumber the satellite galaxies detected in the vicinity of 
the Milky Way and Andr omeda constitutes the so-called 
"miss ing satellite problem" l|Klypin et al ] |l999t lM oore et al.l 
Il999r i. While it has been argued that astrophysical processes 
that quench star for mation in low-mass halos may explai n 
this discrepancy fe.g. iMaccio et alJliuTol : iFont et al.ll201ll ). 



E-mail: ez@astro.su.se 



this implies that a vast population of extremely faint or com- 
pletely dark substructures should be awaiting discovery in 
the halos of galaxies. Provided that CDM is in the form of 
Weakly Interacting Massive Particles (WIMPs), these sub- 
halos are in principle detectable by the Fermi Gamma-ray 
Space Telescope because of their annihilation fluxes. How- 
ever, Fermi has so far fail ed to detect any unambigious sig- 



nal from such objects (e. g. Belikov et al.ll2011 [ Zechlin et al 



20121: |H ooper fc Lindenl 



2012. but see 



Bringmann et al 



2012bl and lSu fc Finkbeineijl2012l for a different view 



Gravitational lensing may provide an independent test 
for the presence of dark ha lo substructures (for a review, see 
IZackrisson fc Riehndl201oT ). A foreground galaxy that hap- 
pens to be aligned with a background light source can pro- 
duce multiple images of the background object, with a typi- 
cal image separation of ~ 1" (an effect known as strong lens- 
ing or macrolensing) . While simple, smooth models of galaxy 
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lenses are usually able to reproduce the positions of these 
macroimages, their observed flux ratios are more difficult to 
explain. Such flux-ratio violations have been interpreted as 
evidence of substantial small-scale structure w i thin the main 
lens ( e.g. iMao fc Schneiderll 19981 ; IChiball20oj iKeeton et ail 
l2003l : lKochanek fc Dara]H2004T ). A notable problem with this 
picture is that current CDM simulations predict too little 
subst ructure to explain many of these flux-rat i o violations 
(e.g.. LMaccio fc Mirandal l200rj: I Xu et al 1 120091 : IChen et all 
2011 - but see iMetcalf fc Amarall2010l ~ possibly pointing to 
a considerable contribu tion from low-m ass halos elsewhere 
along the line of sig ht (|Xu et al.ll2O10h or some additional 
form of substructure (dark or luminous) within the lens. 



A slightly different lensing approach exploits the small- 
scale distortions that halo substructure is expected to in- 
troduce in the morphologies of extended macroimages. Sub- 
structures of mass > 10 s Mq can perturb gravitational arcs 
and Einste in rings on scales resolvable w ith the Hubble Space 
Telescope (|Vegetti fc Koopmandl2009al lbT) and detections of 
~ 10 9 -10 10 Mq obje cts have already been made this way 
l|Vegetti et all l2010al lbl. 120121 ). In line with the flux ratio 
anomaly results, these observations seem to suggest a sub- 
halo mass fraction that is significantly higher than predicted 
by standard C DM, and possibly als o a fla tter subhalo mass 
function slope jVegetti et al.ll2010bl . [ioT^ V 



By mapping extended macrolensed sources with mil- 
liarcsecond or sub-milliarcsecond resolution using Very Long 
Baseline Interferometry (VLBI) techniques at radio wave- 
lengths, substructures at even lower masses can in principle 
be detected. Such objects may i ntroduce kinks and bends 
in multiply-imaged quasar je ts l|Wambsganss fc Paczvnskl 
ll992l ; lMetcalf fc Madaull200ll ) and one detection of a ~ 10 - 
10 7 Mq object has already been claimed using this technique 
iMetcalJ l2002h . In this situation, the lensing effects pro- 
duced by halo substructures can be separated from intrinsic 
morphological features in jets, since the latter would be re- 
produced in all macroimages whereas dark matter clumps 
in the halo of the lens would affect each macroimage dif- 
ferently. Similar methods for exploiting the lensing effects 
produced by halo substructures on milli- and submilliarc- 
second scales have also been explo re d by lYonehara et al.l 
ifeOO^Ifrioue fc Chibal(|2003l.l2005alfB) ; lHisano et all l|200fj) ; 
lOhashi et alj (|2009I ) and iRiehm et all (|2009T l. However, ef- 
fects of this type tend to be sensitive to the density pro- 
files of substructures, and may be undetectable for all 
but the very densest, most extreme forms of substructure 
i|Zackrisson et al.ll2008h . 



Here, we use lensing simulations to explore the 
prospects of using macrolensed quasar jets observed at sub- 
milliarcsecond resolution, in searches for standard CDM 
subhalos, ultracompact minihalos and primordial black holes 
within the main lens. These different forms of substructure 
are described, along with previous constraints on such ob- 
jects, in Sect. [2] The details of our simulations and assump- 
tions are covered in Sect. [3] In Sect, f?] we present our re- 
sults and in Sect.[5]we discuss some lingering issues with the 
adopted technique. Sect. [6] summarizes our findings. 



2 DIFFERENT FORMS OF HALO 
SUBSTRUCTURE 

2.1 Standard CDM subhalos 

At z = 0, the CDM scenario predicts the existence of dark 
matter halos with masses ranging from ~ 10 15 Mq down to 
the cutoff in the density fluctuation spectrum, which is set 
by the detailed properties of the CDM particles. For many 
types of WIMPs, t his cut-off lies so mewhere in the range 
~ 1O _11 -1O~ 3 M (|Bringmannll2009l ). However, alternative 
models involving superweakly interacting particles (super- 
WIMPS) and MeV mass dark matter may p l ace the cutoff 
as hi gh as ~ 10 3 -10 8 M© l|Hisano et al.ll2006l ; iHooper et ail 
l2007f ). 

As these low-mass halos merge to form more massive 
ones, some temporarily survive in the form of subhalos 
within the larger halos. N-body simulations indicate that 
the subhalos within a galaxy-sized CDM halo should follow 
a mass function of the type: 



dN 
dMsub 



oc M~ 



(1) 



with ct ~ 1.9 |Springel et alJl200St iGao et all lioITh . The 
relative contribution from subhalos with mass M > 10 5 Mq 
to the dark matter surface mass density at the typical loca- 
tion of macroim ages in a galaxy-mass lens is / BU b ~ 0.002 
llXu et all [20101 ). albeit with a large scatter (|Chen et al.l 
12QU!). 

The density profiles of isolated field halos in CDM sim- 
ulations c an be reasonably well described b y Navarro, Frenk 
& White H avarro, Frenk. fc White! Il996l ; hereafter NFW) 
profiles: 



P {R) 



Pi 



(R/Rs)(l + R/Rs) 



(2) 



where Rs is the characteristic scale radius of the halo. The 
slope of the inner density cusp (/? = dlnp/dlnr) in this 
profile is f3 — —1, and this makes it difficult for NFW ha- 
los in the dwarf-galaxy mass range to produce millilensing 
effects of the type we are considering in this paper. Typ- 
ically, cusp slopes obe ying /? < —1.5 would be required 
|Zackrisson et al1l2008l ). Later work has shown that models 
with inner cusp slopes that become progressively shallower 
towards the centre provide even better fits to isolated ha- 
los in CDM simulations, eventua lly reaching inner slopes of 
P ^ — 1 (e.g. iNavarro et al-lEjOlCf ). In the context of detect- 
ing millilensing effects from low-mass halos, this just makes 
matters worse, since the central density is reduced. 

Once a halo falls into the potential well of a larger halo 
and becomes a subhalo, it is stripped of material - primar- 
ily from its outskirts - due to interactions with its host halo 
and with other subhalos. This alters the density profile of the 
subhalo compared to an isolated halo of the same mass (e.g. 
lHavashi et al.l2 003: Kaz antzidis et ail l2004), but these mod- 
ifications tend to diminish rather than enhance the ability 
of a CDM subhalo to p roduce detectable millilensing effects 
(|Zackrisson et al.ll2008l '). 

To demonstrate that standard CDM subhalos do not 
provide a significant "background" of millilensing events in 
the observational situation that we consider, we have ther- 
fore adopted NFW profiles for these objects, since this re- 
sults in an overoptimistic estimate on the millilensing effects 
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that standard CDM subhalos are likely to produce. Even 
then, the chances of detecting millilensing effects from these 
objects turn out to be negligible in the observational situa- 
tions we are considering. 

To derive the Rs values of our NFW subhalo profiles, 
we adopt the mass-d ependent concentratio n p arameters c = 
R viT /R s from either iBullock et ail |200ll ) or iMaccio et ail 
(2008), where i? v h- is the virial radius of the halo. Since both 
of these recipes predict higher concentration parameters for 
low-mass halos, and since more centrally concentrated pro- 
files (i.e. profiles with higher c) are more efficient in pro- 
ducing millilensing effects, we calculate the subhalo con- 
centration parameters based on their current masses rather 
than the masses they had prior to becoming subhalos. Since 
nearly all subhalos have lost co nsiderable amounts of ma- 
terial fe.g. lVale fc Ostrikedliooll . this also results in overly 
optimistic millilensing properties. 



2.2 Intermediate-mass black holes 

Intermediate-mass black holes (IMBHs; here assumed to 
have masses ~ 10 3 -10 6 Mq) may either form primordially 
(typically when the Universe is <C 1 s old), or due to the 
collapse of baryonic objects later on. The primordial variety 
could in principle comprise a substantial fraction of the dark 
matter, although a host of observational constr aints makes 
this s eem unlikely (for a recent compilation, see lCarr et al.l 

I2010I ). 

The strongest upper limits on the cosmological den- 
sity of primordial black holes in the 10 3 -10 5 Mq mass 
range come from the effect that accretion onto these objects 
would have on the cosmic micro wave background radiation 
(|Ricotti. Ostriker fc Mack! 120081 ). Primordial black holes 
with masses ~ 10 3 -10 4 Mq are also strongly constrained by 
the effects of gravity waves on pulsar timin g measurements 
l|Saito fc Yokovamai |2010| ; ICarr et al.l |2010|). and at M > 
10 4 Mq by dynamical constraints l|Carr fc Sakcllariadou 



Kassiola et al 



1999). Using a technique first proposed by _ 
l|l99ll ). IWilkinson et all (|200lT ) moreover used the absence 
of millilensing effects in non-macrolensed radio sources to 
place upper limits on IMBHs at M > 10 Mq. It has, how- 
ever, been argued that some of these constraints may be 
sidestepped under certain circumsta nces, and that bot h the 
size evolution of early-ty pe galaxies jTotani et al.ll2010l ) and 
entropy considerations (jFrampton et al.l |2010| ) favour sce- 
narios in which essentially all of the dark matter is in the 
form of ~ 10 Mq primordial black holes. 

Intermediate-mass black holes that were not produced 
primordially may instead either form as the end products 
of ve ry massive population III stars (e.g. iMacieiewski et all 
l200ll ). through the direct collapse of gas in small halos 
at high redshift (e.g. iBegelman et al l T2OO6I) or the col- 
lapse of dense star clusters fe.g. lDevecchi fc Volonteri 2009; 
iDavies et al.ll201ll ) . Such IMBHs may no w be hiding in glob- 
ular clusters_(e ; g_ ; _^^s£erin^^^] [20l3), in satellite galax- 
ies (|van Wassenhove et al l |201oFT or be fre ely floating in 
the halos of galaxies (e.g. iMicic et al.ll201ll ). There is in- 
deed some eviden ce for IMBHs in globular clusters (e.g. 
iNovolaet ai]|2010l ). and IMBHs may also explain some of 
the ultr aluminous X-ray sources detected in o ther galax- 
ies (e.g. iFeng fc Sorial l201ll ; IWebb et all |2012P . However, 
since only a small fraction of the cosmic baryons can be 



locked up in these non-primordial IMBHs, their relative con- 
tributions to the halo masses of galaxies are typically ex- 
pected to be small (/i mbh < 10~ 5 ; e.g. Ilslam et al.l 12004 ; 
iKawaguchi et al.ll2008l ). 

When simulating the potential millilensing effects of 
IMBHs, we treat the surface mass density fraction /imbh 
in IMBHs at the position of the macroimages as a free pa- 
rameter, and for simplicity assume all IMBHs to have the 
same mass. In the case where the number density profile 
of IMBHs has the same shape as the density profile of the 
dark halo, /imbh also corresponds to the halo mass fraction 
in IMBHs. 

In principle, primordial black holes may over time ac- 
crete substantial amounts of dark matter and develop dark 
matter halos of their own l|Mack et al.l 120071 1. similar to 
the ultracompact minihalos discussed in Sect. 12.31 IMBHs 
forming through the collapse of pop III stars in miniha- 
los may also be su rrounded by their ow n, highly contracted 
dark matter halos |Sandick et al.|[201ll ) . Such compound ob- 
jects are expected to have lensing properties intermediate 
between IMBHs and ultracompact minihalos, but are not 
treated in detail in our simulations. 



2.3 Ultracompact minihalos 

Primordial density perturbations with Ap/p = 8 < 0.3 are 
too small to form primordial black holes as they enter the 
horizon. Those with 8 > 10 -3 may nonetheless still be large 
enough to seed the formation of ultracompact minihalos 
(UCMHs). Such perturbations might be produced during 
phase transitions, around topological defects, or in the pri- 
mordial spectrum of perturbations from inflation. The dark 
matter contained in these perturbations would collapse into 
UCMHs shortly after matter-radiation equality, via radial 
infall from a universally cold, smooth cosmological back- 
ground. This radial collapse would leave UCMHs with much 
steeper central density profiles than standard CDM halos 
jRicotti fc Gouldll2009l ). 

If dark matter exists in the form of self-annihilating 
WIMPs, UCMHs would be gamma-ray emitters, and strong 
limits on their cosmological density have already been 
derived from the effect that this would have on Fermi- 
LAT source identifications, the diffu se gamma-ray back- 
ground and co s mic r ei onization dScott fc Sivertssod 120091: 
Josan fc Greenl |20ld; lYang et all l2011al lbl: IZhand l201ll ; 
Bringmann et alll2012a ). If dark matter does not annihilate, 



UCMHs in the ~10 -10 Mq range may still be detectable 
in the futu re by their as trometric lensing effects on Milky 
Way stars (|Li et al.1120121 ). 

Here, we explore to what extent submilliarcsecond 
observations of macrolensed jets would be able to con- 
strain the properties of UCMHs. Because dark matter self- 
annih ilation would reduce th e central density of UCMHs 
(e.g. IScott fc Sivertsson|[2009l ). UCMHs made out of self- 
annihilating WIMPs would not be efficient millilenses. We 
therefore focus on UCMHs made out of non-annihilating 
dark matter (e.g. asymmetric dark matter, axions, sterile 
neutrinos). 

Radial infall leads to a density profile p cx r~ 2 ' 25 , 
slightly steeper than the par -2 profile of a singular isother- 
mal sphere (often used to model lensing by the inner regions 
of large galaxies). The dark matter profile in a UCMH (see 
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iRicotti fc Gouldl 120091 ; iBringmann et all l2012al , for a de 
tailed discussion) is given by 

3/cdmMucmh(z) 



p(r,z) 



(3) 



167r7?uCMH(z) 4 r* 

where /cdm is the cosmological fraction of matter in CDM, 
Mucmh(z) indicates the UCMH mass at redshift z, and 



-Rucmh(z) 
pc 



0.019 



1000 
z+1 



( Mucmh(z) 



(4) 



is the UCMH radius, defined as the distance within which 
the density is at least twice that of the cosmological back- 
ground. Following matter- radiation equality at z cq , a UCMH 
born from an initial dark matter overdensity of mass Mi ac- 
cretes both dark and baryonic matter as 



Mucmh(z) = ^^±±Mi 

2+1 



(5) 



This accretion presumably cuts out when the cosmological 
background is no longer smooth, i.e. when the first substan- 
tial structure formation occurs and the smallest star-forming 
minihalos appear. In this case, present-day UCMH masses 
and radii can be obtained by setting z ~ 10 in Eqs. ([3| and 
©, so that 

mS C mh = MucmhO < 10) = MucmhO = 10), (6) 
i?ScMH = Rvcmh(z < 10) = JJtjcmhO = 10). (7) 

However, the finite temperature of the smooth cosmo- 
logical background from which UCMHs accrete softens the 
density profile in the innermost region, due to conservation 
of angular momentum. This can be conservatively modelled 
as a cutoff at some inner radius r m i n , inside which one as- 
sumes the density to be constant. Following previous work 
(Brin gmann et al. 2012a), we adopt this strategy here, tak- 
ing a flat density profile within 



-^iicmh 



2.9 x 10" 7 



1000 

z c + 1 



i "UCMH 



(8) 



Here z c refers to the redshift of UCMH collapse (the point 
at which the growth of the matter overdensity becomes non- 
linear); we adopt z c = 1000, also in line with e arlier work 
ijRicotti fc Gouldll2009l : IBringmann et al.ll2012ah . 



3 LENSING SIMULATIONS 

To simulate the effects of dark halo substructure on the mor- 
phologies of macrolensed jet s, we use a numer i cal sc heme 
similar to that developed by iMetcalf fc Madaul (|2001T) . An 
extended source is assumed to be multiply-imaged by a fore- 
ground galaxy, and the lens equation is used to determine 
the lens plane positions of the corresponding macroimages. 
A small region around each such macroimage is then popu- 
lated with randomly distributed dark halo substructures and 
simulated in greater detail. The deflection angles (with con- 
tributions both from substructures and the macrolens) are 
computed for every pixel within this region and converted 
into a numerical surface brightness map of the macroimage. 
These maps are initially generated with a very fine pixel 
scale, but are then convolved with a Gaussian filter to match 
the finite resolution of the VLBI arrays we consider. Both 



the resolution and the intrinsic source dimensions are deter- 
mined by the frequency at which we assume the jets to be 
observed, as described in Sect. I3.ll and l3.2l 

The macrolens is modelled as a singular isothermal 
sphere (as appro priate for early-type galaxies a cting as 
strong lenses; e.g. iRusin. Kochanek fc Keetonll2003h at Z\ = 
0.5, with line-of-sight velocity dispersion a v — 240 km s _1 , 
giving a mass of ~ 10 13 Mq within the virial radius, and 
two macroimages with separation w 2". We furthermore ad- 
just the alignment of the source and main lens to ensure 
a macrolens magnification that is not unrealistically high. 
The simulations presented in this paper are all based on a 
lens-source configuration that in the absence of substructure 
would give magnifications fii ~ 10 and £t2 ~ 8 for the two 
macroimages. All simulations are based on a ACDM cosmol- 
ogy with H = 70 km s" 1 Mpc" 1 , fi M = 0.27, Q A = 0.73. 

When distributing halo substructures within the simu- 
lated region, we for simplicity assume that the surface mass 
density across the macroimage is completely dominated by 
dark matter. While this assumption may be violated in 
multiply-imaged systems where one of the macroimages hap- 
pens to lie very close to the lensing galaxy, this nonetheless 
seem s to be a fair approx i mation in the maj ority of cases 
(e.g. iBegelman et al.fl201ll ; iPoolev et al]|2012l ). In the case 
of intermediate mass black holes and ultracompact miniha- 
los, we moreover assume that their number densities trace 
that of the dark matter. The surface number density of such 
subtructures then simply depends on their relative contri- 
bution to dark matter f2 S ub/f2cDM an d their mass distri- 
bution. These dark matter fractions in intermediate-mass 
black holes and ultracompact minihalos are referred to as 
/imbh and /ucmh respectively. Since detailed predictions 
for the mass distribution of IMBHs and UCMHs are highly 
model-dependent, we assume all such objects to have the 
same mass, which we then vary to explore what parts of the 
(/imbh,Mimbh) or (/ucmh,Mucmh) parameter space that 
a given set of observations would be able to probe. 

In the case of standard CDM subhalos, we adopt the 
mass distributions inferred from either simulations or obser- 
vations. As discussed in Sect. 12.11 current simulations sug- 
gest /nfw = 0.002 at the typical positions of macroimages 
in galaxy-sized dark halos, and the subhalo mass function 
given by Eq. {T}. H owever, since the recent len sing detec- 
tions of subhalos by IVegetti et ail l|2010bl . 120121 ') hint at a 
flatter mass function and a mass fraction that is an order of 
magnitude higher, we also explore the consequences of set- 
ting /nfw = 0.03 and changing the mass function slope of 
Eq. Qtoa = 1.1. 



3.1 VLBI observations 

A number of macrolensed radio jets are already known 
and have been observed using VLB I technique s (e.g . 



Rusin et all |2002| ; iBiggs et al.l 12004 lYork et all 120051 ). al- 



mg 

Garrett et all Il994 [King et all Il997l; IRos et all 



2000; 



though typically not with arrays capable of resolving sub- 
milliarcsecond scales. Designing a survey aimed to search for 
small-scale distortions in targets like these does, however, 
also involve other considerations than just the resolution. 
The frequency at which one chooses to observe these jets 
limits the resolution at which the jets can be mapped using 
suitable VLBI arrays, but also affects the intrinsic source 
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Figure 1. Simulated radio maps of strongly lensed quasar jets at 86 GHz, 22 GHz and 8.4 GHz, respectively from left to right (source 
sizes 2 X 0.5 pc, 10 X 2.5 pc and 40 X 10 pc), subject to macrolensing by the main lens. The two subplots of each image show the 
two macroimages of the source. The bottom row contains the contour representations of the macroimages in the upper row, with the 
outermost contours corresponding to 10% of the peak flux in these images. Please note the different scales of the images at the three 
frequencies. 



size (e.g. iTorniainen et al1l2008h . To identify the observa- 
tional strategy that maximizes the scientific output in terms 
of detection prospects for dark halo substructure, we here 
consider observations at three different frequencies, each us- 
ing a different VLBI array: 

• Observations at 8.4 GHz using theglobal array, includ- 
ing the European VLBI Network (EVrQ) and the Very Long 
Baseline Array (VLBA0), giving a resolution of ~ 0.7 mil- 
liarcseconds 

• Observations at 22 GHz using the EVN, giving a reso- 
lution of w 0.25 milliarcseconds 

• Observations at 86 GHz using the full Atacama Large 
Millimeter (ALMAn) array (66 antennas) connected to the 
global 3 mm arrajQ giving a resolution of « 0.05 milliarcsec. 
This observing mode is not available at the current time, but 
is likely to come on line in a few years. 

These arrays also have different sensitivities, which con- 
strains the numbers of potential targets and also the ap- 
parent lengths of the jets. However, since we are simulating 
the effects of generic sources rather than individual targets, 
this is not addressed in our current simulations. 



1 http://www.evlbi.org/ 

2 http://www.vlba.nrao.edu/ 

3 http:/ /www. almaobservatory.org/ 

4 http:/ /www. mpifr-bonn.mpg.de/div/vlbi/globalmm/ 



3.2 Source size and morphology 

We assume the source to be an intrinsically straight jet 
with a 2-dimensional Gaussian surface brightness profile and 
length 40, 10 and 2 pc at 8.4, 22 and 86 GHz, respectively, 
and a width that is a quarter of the length. These sizes are in 
a r easonable agreem e nt wit h source size estimates presented 
in ITorniainen et all l|200Sf ) and the jet lengt hs calculated 
from the MO J AVE sample l|Lister et al.ll2009l fl. 

The source morphology and length-to-width ratio are 
mainly adopted for illustrative purposes. The limits pre- 
sented in Sect. [4] are fairly insensitive to these assumptions, 
and depend mainly on the intrinsic jet area. At fixed angu- 
lar resolution, a larger jet results in a stronger constraint 
whereas a smaller jet makes the constraints weaker. The 
results of Sect. [4] can therefore be rescaled to accommo- 
date other jet sizes. In broad terms, our assumptions on 
the so urce sizes are similar to those used by llnoue fc Chibal 
(|2003h . 



4 RESULTS 

In Fig. [T] we present our simulated images of strongly lensed 
quasar jets at 86, 22 and 8.4 GHz. Each image pair in the fig- 
ure corresponds to the two macroimages of a single radio jet 
as produced by the main lens in the absence of any millilens- 
ing effects. While initially generated using a much smaller 
pixel scale, these images have been degraded using a Gaus- 
sian filter to match the resolution relevant for observations 

5 http://www.physics.purdue.edu/MOJAVE/ 
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Figure 2. Examples of simulated radio maps of a macrolensed quasar jet at 86 GHz (source size 2 X 0.5 pc and resolution 0.05 
milliarcseconds), subject to millilensing distortions by IMBHs with /imbh = 0.02 in the halo of the main lens. Each image pair represents 
the two macroimages from Fig.[T] distorted by millilensing effects from IMBHs with either Mjmbh = 10 3 i 10 4 j 10 5 ^© or 10 6 Mq. The 
positions of the IMBHs are indicated by red dots. The slight macroimage distortions and displacements seen in right panels for the 10 5 
or 10 6 Mq cases are produced by IMBHs just outside the plotted region. 




mas mas mas mas mas mas 

Figure 3. Examples of simulated radio maps of a macrolensed quasar jet at 22 GHz (source size 10 X 2.5 pc and resolution 0.25 
milliarcseconds), subject to millilensing distortions by IMBHs with /imbh = 0-01 in the halo of the main lens. Each image pair represents 
the two macroimages from Fig.[T] distorted by millilensing effects from IMBHs with either Mjmbh = 10 4 , 10 5 or 10 6 Mq. The positions 
of the IMBHs are indicated by red dots. 

at these frequencies (0.05, 0.25 and 0.7 milliarcseconds re- ures depict how these contours are distorted by various kinds 
spectively). The bottom row shows the corresponding isofiux of halo substructure within the main lens, 
contour plots, where the outermost contours correspond to 
w 10% of the peak flux in these images. All subsequent fig- 
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Figure 4. Examples of simulated radio maps of a macrolensed quasar jet at 8.4 GHz (source size 40 X 10 pc and resolution 0.7 
milliarcseconds), subject to millilensing distortions by IMBHs with /imbh = 0.005 in the halo of the main lens. Each image pair 
represents the two macroimages from Fig. [T] distorted by millilensing effects from IMBHs with either Mimbh = 10 s or 10 6 Mg. The 
positions of the IMBHs are indicated by red dots. 



4.1 Detecting intermediate- mass black holes 

In Fig. [2] we present examples of the simulated macroimages 
in the case where a fraction /imbh = 0.02 of the dark halo 
of the main lens is in the form of intermediate-mass black 
holes with mass Mimbh = 10 3 — 10 6 Mq. In this case, 86 
GHz observations (ALMA + global array) are assumed, im- 
plying the smallest jet size (intrinsic length 2 pc) and the 
highest resolution (0.05 milliarcseconds) considered in this 
paper. Since /imbh is kept fixed, the number of IMBHs per 
unit area drops by a factor of 10 3 in the lens plane when 
going from Mimbh = 10 3 Mq to 10 6 Mq. However, because 
the more massive IMBHs also have larger Einstein radii, po- 
tentially detectable distortions are produced in all the cases 
plotted. Since the distortions in the two macroimages are 
uncorrelated, millilensing should also be straightforward to 
separate from intrinsic jet features (but see Sect. [S] for po- 
tential caveats). 

The probability of seeing millilensing effects in at least 
one macroimage of a given two-image system depends on 
Mimbh, the angular resolution and /imbh, but is deemed to 
be P m iiii > 50% in all the simulations presented in Fig. [2] 
When analysing a survey of N such macrolens systems, 
the probability Pdetection of detecting millilensing becomes 
^detection = 1 - (1 - P m mi) N - By adopting P mim > 50%, one 
should therefore be able to rule out /imbh = 0.02 for IMBHs 
in the mass range Mimbh = 10 3 -10 6 M at the > 95% level 
by surveying N ~ 5 systems. By further increasing the size 
of the survey, even lower /imbh can in principle be probed. 
To first order P m mi scales with /imbh, so that P m mi > 5% 
if /imbh ~ 0.002. Hence, to reach a detection probability 
of Pdetection > 68% if /imbh = 0.002, one needs to observe 
N > 22 systems. 

As demonstrated in Fig. [3] the larger source size as- 
sumed for the 22 GHz (EVN) observations (10 pc) allows 
IMBHs with dark matter fractions as low as /imbh = 0.01 
to be detected with probability P m mi > 50%. The lower 
resolution (0.25 milliarcseconds) provided by the EVN at 
the same time prohibits the detection of IMBHs with mass 
Mimbh ~ 10 3 Mq. By surveying N ~ 5 systems, one 



should be able to rule out /imbh = 0.01 for Mimbh ~ 10 4 - 
10 6 Mq at 95% confidence level. A detection probability of 
Pdetection > 68% can also be reached at /imbh = 0.001 if one 
is able to observe N ^ 22 systems. 

Similarly, the even larger jet (intrinsic length 40 pc) 
adopted for our simulated 8.4 GHz observations allows for 
stronger constraints on /imbh, but the lower resolution (0.7 
milliarcseconds) at the same time limits the IMBH mass 
range for which millilensing effects can be detected. Still the 
macroimage distortions produced by 10 5 -10 6 Mq IMBHs 
would be detectable with this resolution, and such effects 
would turn up with P m mi > 50% probability even if the 
IMBH halo mass fraction is as low as /imbh = 0.005. Fig. [4] 
includes an example of such millilensing distortions pro- 
duced by 1O 5 M and 10 6 M Q IMBHs with /imbh = 0.005. 

Table [1] summarizes the /imbh limits that observations 
of a single macrolensed jet (one image pair) at 86, 22 and 
8.4 GHz would be able to probe (with ^ 50% detection 
probability). As described above, these limits can easily be 
rescaled to accommodate observations of larger number of 
multiply-imaged systems. The constraints resulting from a 
survey of N ~ 5 systems would p roduce constraints that 
are a factor of a few better than the lWilkinson et al.l (2001) 
millilensing constraints on 10 6 Mq primordial black holeqj. 
There are no competitive lensing constraints on IMBH at 
10 3 -10 5 Mq, but there are still a host of other constraints 
that may be applicable, in particular those r elated to accre- 
tion onto these objects fsee lCarr et al.l 2010. for a review). 



6 Formally, the Wilkinso n et all l l200ll ) constraints apply to 
millilcnses located anywhere along the line of sight to the ra- 
dio sources in their sample (mean redshift z as 1.3), whereas ours 
apply only to millilenses within the main lens. The difference may 
be relevant in scenarios in which the IMBHs do not follow the dis- 
tribution of dark matter on large scales (e.g. if they are formed 
through baryonic processes in the vicinity of galaxies.). Moreover, 
our approach could in principle produce somewhat stronger con- 
straints if we were to consider millilenses along the entire line of 
sight 
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Figure 5. Examples of simulated radio maps of a macrolensed quasar jet at 86 GHz (source size 2 X 0.5 pc and resolution 0.05 
milliarcseconds), by UCMHs with Mtjcmh = 10 6 and 10 s Mq in the halo of the main lens. The probabilities of detecting such effects 
are, however, negligibly small unless the UCMH dark halo fraction is /ucmh ~ 0.2. 



Table 1. The lowest halo mass fraction in IMBHs, /imbHi that 
would produce detectable millilensing distortions with P m jnj > 
50% probability in a single macroimage pair. 



Frequency (GHz) 


Mass (Mq) 


Source size (pc) 


min /imbh 


86 


10 3 - 10 6 


2 X 0.5 


2 X 10~ 2 


22 


10 4 - 10 6 


10 X 2.5 


1 x 10~ 2 


8.4 


10 5 - 10 6 


40 x 10 


5 x 10~ 3 



By probing N ~ 11 (28) macroimage pairs, the detec- 
tion probability can be pushed to Pdetection ~ 68% (95%) 
at these /ucmh limits. In order to probe UCMH dark halo 
fractions significantly below /ucmh ~ 0.1, hundreds of im- 
ages would therefore need to be observed. While there are no 
competit ive lensing const r aints at ~ 10 6 Mq, it is possible 
that the IWilkinson et al.l l|200ll ) observations of 300 z ~ 1 
radio sources (not macrolensed) at ~ 1 milliarcsecond res- 
olution would be able to do better for ~ 10 s Mq UCMHs 
than the predicted limits we give in this paper. 



Table 2. The lowest halo mass fraction in UCMHs, /uCMHi that 
would produce detectable millilensing distortions with P m jui > 
10% probability in a single macroimage pair. 

Frequency (GHz) Mass (Mq) Source size (pc) min /tjcmi 



86 
22 
8.4 



10 6 
10 7 
10* 



10 8 
10 8 



2 x 0.5 
10 x 2.5 
40 x 10 



2 x lO" 1 
1 x lO" 1 
5 x 10~ 2 



4.2 Detecting ultracompact minihalos 

When compared to IMBHs of the same mass, UCMHs have 
much smaller Einstein radii and are far more difficult to 
detect through millilensing effects. Our simulations show, 
that while Mtjcmh ~ 10 6 -10 8 Mq UCMHs may in principle 
be detectable through small-scale macroimage distortions, 
the probability of observing this effect is exceedingly small 
unlessthe UCMH dark matter fraction /ucmh is very large. 

In Fig. [5] we show examples of the millilensing distor- 
tions that 10 6 and 10 s Mq UCMHs would produce in the 
case of 86 GHz observations (ALMA + the global array). 
However, the probability of seeing effects of this type in 
a given macroimage pair is only P m nii ~ 10%, even for a 
UCMH dark halo fraction as high as /ucmh = 0.2. The 
detection prospects become somewhat better at 22 and 8.4 
GHz (/ucmh > 0.05-0.1 at P mim « 10%) due to the larger 
source adopted sizes at these frequencies, but only for 10 7 - 
10 8 M e objects (see Table [2]). 



4.3 Detecting standard CDM subhalos 

As expected, standard CDM subhalos (assumed to have 
NFW density profiles) are not detectable using the obser- 
vational scheme considered in this paper. In Fig. [6] we show 
examples of the ~ 10 7 and ~ 10 8 Mq subhalos close to the 
macroimages i n our simulations at 22 GHz, assuming the 
iBullock et alJ (|200lh relation between mass and concentra- 
tion parameter. In the standard case of / su b ~ 0.002, only 
M < 1O 7 M0 subhalos are sufficiently numerous to have 
a decent probability of showing up in the vicinity of the 
macroimages, and even though such objects may affect the 
overall magnification and the curvature of the jet (see left 
panel of Fig. [6|, the associated small-scale distortion is too 
small to be resolved. NFW subhalos at ~ 10 8 Mq may in 
principle give rise to the detectable distortions (right panel 
of Fig. [6}, but the probability of attaining the required alig- 
ment between macroimage and subhalo is negligibly small. 
We estimate that the probability of detecting small-scale dis- 
tortions due to 1O 8 M0 NFW subhalos in a single macroim- 
age pair is no more than P m im ~ 3 x 10 -4 in this case. At 
10 9 M Q , the probabilit y is even lowe r (P m m i ~ 4 x 10~ 5 ). 



Recent results bv lvegetti et alJ (|2010bl . [2012] ') hint at a 
higher surface mass density contribution (/ su b ~ 0.03) and 
a flatter subhalo mass function (a = 1.1 in Eq. JIJ) than 
predicted by current CDM simulations, but even if we adopt 
these values, the probability for detection remains too low 
(» 5 x 10~ 4 for NFWs of mass 10 8 -10 9 Mq) to make this 
search strategy attractive. 
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Figure 6. Examples of simulated radio maps of a macrolensed quasar jet at 22 GHz (source size 10 X 2.5 pc and resolution 0.25 
milliarcseconds), subject to millilensing distortions by NFW subhalos with masses of either ~ 10 7 or 10 s Mq. The red dots mark the 
positions of the centres of these subhalos. As the left pair of images shows, ~ 10 7 Mq NFWs produce very mild distortions only, whereas 
NFW subhalos of mass ~ 10 8 Mq may produce more significant distortions if they are placed sufficiently close to a macroimage. The 
probability for such superpositions to occur is, however, very small. 



The problem is one of source size - the macrolensed jets 
we consider cover an area in the lens plane that is several 
orders of magnitude too small to intersect such massive sub- 
halos. The intrinsic source size would need to have an area 
~ 10 3 times greater than the largest jets we consider (40 x 10 
pc at 8.4 GHz) to push the detection probabilities into the 
interesting range P m mi > 10%. This essentially requires a 
completely different ki nd of source, like the d usty sub-mm 
galaxies considered by llnoue" fc Chibal i|2005ah . 

These result are admittedly sensitive to the concen- 
tration parameters adopted for the NFW subhalos. In the 
examples abo v e, we have used the c(M v i r ) relation from 
iBullock et al.l (|200ll ). which for objects in the relevant 
mass range (10 7 -10 8 Mq) results in concentration param- 
eter s a factor of ~ 2 h igher than the ones predicted by 
the Macc io~et al.l (120081 ) relation. If we instead adopt the 
iMaccio et alj ( 20081 ) c(M v w) scaling, the detection thresh- 



old shifts upward by an order of magnitude in mass, so that 
image distortions predicted for 10 7 and 10 s Mq NFW in the 
Bullock et al. case (Fig. [6} instead are produced at masses 
of ~ 10 8 and ~ 10 9 Mq. 



5 DISCUSSION 

In previous sections, we have argued that millilensing- 
induced distortions of quasar jets may be distinguished 
from morphological features intrinsic to these sources, since 
the latter would be reproduced in all macroimages whereas 
millilensing should affect each image independently. How- 
ever, this argument comes with a caveat. The time delay 
between the images in quasar-galaxy lenses can be up t o 
a year (for a compilation of time delays, see lOgurll 120071 ). 
which means that intrinsic, transient features in the jet may, 
at any given time, be visible in just one of the images and be 
mistaken for millilensing effects. This is for instance likely 
to be the case in superluminal radio jets, where blobs are 
seen to move ~ 1 milliarcseconds yr _1 along the jet (e.g. 
Ijorstad et al.ll2~00ll ). For macrolensed jets that show signs 



of millilensing distortions, it may therefore become neces- 
sary to obtain data at two or more epochs. Since halo sub- 
structures give rise to millilensing magnification pattern that 
will appear stationary over decadeqj, any distortions that 
seem to move along the jet are bound to be intrinsic to the 
source. Small-scale features that are not duplicated in the 
other macroimages and appear with a fixed angular position 
(as, for instance, measured from the base of the jet) over the 
course of more than a year is on the other hand likely caused 
by millilensing. 

The detection of milliarcsecond or submilliarcsecond- 
scale image distortions would prove the existence of sub- 
structures within the macrolens, and also allow constraints 
on their surface number densities to be set. However, the 
exact nature of t he millilenses may still be very difficult to 
determine. While llnoue fc Chibal l|2005b| ) has demonstrated 
that the distortions induced in extended images (like the 
ones we model here) contain some information about the 
density profiles of the lenses, the finite resolution and sensi- 
tivity of actual observations could still allow for considerable 
degeneracies in cases where neither the masses nor the den- 
sity profiles of the millilenses are known. In our simulations, 
IMBHs and UCMHs can for instance sometimes produce 
very similar lensing distortions (although at different masses 
- an UCMH typically needs to be a factor of ~ 10 3 more 
massive than an IMBH to reproduce a given feature). 



6 SUMMARY 

Using simulations of strongly lensed quasar jets, we argue 
that very dense forms of halo substructure (intermediate- 
mass black holes and ultracompact minihalos) within the 



7 assuming a typical effective transverse velocit y due to the bulk 
motions of the source, macrolens and observer ( Gil-Merin o et all 
2005), and that the velocity of the millilenses within the 
macrolens is similar to the velocity dispersion of the latter (240 
km s — 1 in our case) 
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main lens may reveal itself through small-scale morpholog- 
ical distortions in the macroimages. Such distortions can 
be distinguished from intrinsic source features by obtain- 
ing data at multiple epochs. By mapping a handful of 
macrolensed jet systems at submilliarcsecond resolution, we 
argue that ~ 10 3 -10 6 Mq intermediate-mass black holes 
can be detected or ruled out if they contribute a surface 
mass fraction of /imbh > 0.01 (depending on the VLBI ar- 
ray and frequency used) to the dark matter of the main 
lens at the macroimage positions. Ultracompact miniha- 
los in the ~ 10 6 -10 8 Mq mass range may similarity pro- 
duce detectable small-scale effects if such objects comprise 
/ucmh > 0.1 of the dark matter. While standard CDM sub- 
halos at masses of > 10 8 M© can in principle also produce 
milliarcsecond-scale distortions, provided that such objects 
are projected sufficiently close to the macroimages, the prob- 
ability of this is too small (P m mi ~ 10 -4 ) for sources of the 
type we consider. 
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